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ABSTRACT 

We present multifiber echelle radial velocity results for 55 1 stars in the Sextans dwarf spheroidal galaxy and 
identify 294 stars as probable Sextans members. The projected velocity dispersion profile of the binned data 
remains flat to a maximum angular radius of 30'. We introduce a nonparametric technique for estimating the 
projected velocity dispersion surface, and use this to search for kinematic substructure. Our data do not confirm 
previous reports of a kinematically distinct stellar population at the Sextans center. Instead we detect a region 
near the Sextans core radius that is kinematically colder than the overall Sextans sample with 95% confidence. 

Subject headings: galaxies: dwarf — galaxies: individual (Sextans) — galaxies: kinematics and dynamics — 
(galaxies:) Local Group 
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1. INTRODUCTION. OBSERVATIONS AND DATA 

The most recent kinematic studies of dwarf spheroidal 
(dSph) galaxies show that radial velocity (RV) dispersion pro- 
files of these dark matter dominated systems are generally 
flat, with some evidence for a declining dispersion at large ra- 
dius in Draco and Ursa Minor (Wilkinson et al. 2004, but see 
Munoz et al. 2005; Walker et al. 2006; Mateo et al. in prep.). 
However, localized regions having enhanced stellar density 
and/or comparatively "cold" kinematics have been detected 
in some dSphs (Olszewski & Aaronson 1985; Kleyna et al. 
2003; Kleyna et al. 2004, "K04" hereafter; Coleman et al. 
2004). Such substructure may be related to the star formation 
history of individual systems (Tolstoy et al. 2004; Olszewski 
et al. 2006), or may hint at a merger history (Coleman et al. 
2004; Wilkinson et al. 2005). The presence of substructure 
in these tiny galaxies has cosmological implications as well. 
Kleyna et al. (2003) demonstrate that the apparent substruc- 
ture in Ursa Minor should quickly disperse in the presence of 
a centrally cusped dark matter halo, but can remain intact for 
a Hubble time in the potential of a cored halo. This raises 
doubts about whether the smallest galaxies conform to the 
universal density profile that results from simulations of cold 
dark matter (Navarro, Frenk & White, 1997). 

Here we present first results from an ongoing survey of 
dSph stellar RVs. We used the Michigan-MIKE Fiber System 
(MMFS) at the 6.5m Magellan/Clay telescope at Las Cam- 
panas Observatory, in March 2004 and February 2005 to ob- 
tain spectra of target stars selected from V,I photometry of the 
Sextans red giant branch. A complete description of MMFS is 
forthcoming (Walker et al. in prep.). Briefly, 256 fibers placed 
over a field of diameter 30' feed the two channels of the MIKE 
echelle spectrograph (Bernstein et al. 2003). Order-blocking 
filters isolate the spectral region 5 140-5 1 80 A, which contains 
strong Mg-I triplet absorption features and is free of the sky 
emission lines that contaminate the infrared Ca triplet. After 
binning, the spectra have resolution R ~ 20000 (blue channel) 
and R ~ 15000 (red channel). 
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FIG. 1 . — left: Heliocentric RV vs. angular distance from the Sextans center 
(«2000 =10:13:03, <5 20 oo =-01:36:54; Mateo 1998). Dotted, dashed, and long- 
dash lines mark boundaries of N=276, N=294, and N=303 member samples, 
respectively. Not shown are 21 observed stars at 60' < R < 80', none of 
which are likely members, right: Histogram of the Sextans RV distribution. 



The data reduction procedure is similar to that described in 
Mateo et al. (2006, in prep.). Briefly, after extracting and cali- 
brating spectra, we measure RV using the IRAF task FXCOR, 
cross-correlating the spectrum of each target star against that 
of a high-S/N template composed of co-added spectra of RV 
standard stars. A series of quality controls yields a data set 
containing 810 RV measurements for 551 Sextans candidate 
members. We use multiple measurements of 173 stars from 
12 partially-overlapping fields to estimate measurement un- 
certainties as in Walker et al. (2006). The mean la RV uncer- 
tainty is ±2.6 km s" 1 . 

Figure^plots the measured heliocentric RV versus angular 
distance from the Sextans center. The broad distribution of 
RV outliers away from the main peak at 226 km s" 1 is consis- 
tent with the RV distribution of foreground Milky Way dwarfs 
(whose magnitudes and colors place them inside our Sextans 
RGB selection region) predicted by the Besancon Milky Way 
model (Robin et al. 2003) 4 . To determine membership we use 
a biweight estimator (Beers et al. 1990) of distribution cen- 
ter and variance, and iteratively reject RVs that lie more than 
3cr from the distribution center. This algorithm yields a de- 
fault sample of 294 probable Sextans members. Because the 
conventional 3cr cutoff is somewhat arbitrary, we repeat our 
analyses for possible samples containing 276 stars (2.6a out- 
lier rejection), and 303 stars (4a). 



see http://bison.obs-besancon.fr/modele 
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FIG. 2. — Sky map, in standard coordinates, of 303 probable Sextans mem- 
ber stars with velocities measured using MMFS. Color scale indicates the 
magnitude of each stellar RV relative to the sample mean. The ellipse has ge- 
ometric mean radius corresponding to r c ore = 16.6' ± 1 .2' (Irwin & Hatzidim- 
itriou 1995). The dashed (dotted, dot-dash) circle encloses the region from 
the N=294 (N=276, N=303) sample where the measured RVs are most likely 
to have been drawn from two distinct Gaussian distributions. 
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FIG. 3. — Sextans projected velocity dispersion profile measured from the 
N=276, N=294, and N=303 member samples. Bins in (a) have regularly- 
spaced outer boundaries. In b-c, bins contain approximately equal numbers 
of stars. Horizontal error bars indicate the standard deviation of R values. 



2. RADIAL VELOCITY DISPERSION: PROFILE AND SURFACE 

Figure[2]maps the projected positions of the likely Sextans 
members and indicates the magnitude of each star's RV rel- 
ative to the sample mean. We detect no rotation — the max- 
imum rotation signal of 0.9 km s" 1 is exceeded in 43% of 
Monte Carlo trials. Following the procedure described in 
Walker et al. (2006), we estimate from the default N=294 
sample a mean RV of 225.8 ±0.5 km s" 1 and RV dispersion 
8.9 ±0.4 km s" 1 , in agreement with K04. For N=276 and 
N=303 samples the RV dispersion is 7.1 ±0.3 and 10.3 ±0.5 
km s , respectively. Binning by radius, we estimate the 
projected RV dispersion profile, <r p (R), using the Gaussian 
maximum-likelihood method of K04 (see their Equation 1). 
The RV dispersion calculated for a given bin is not affected 
if we modify K04's Equation 1 to include a fractional com- 
ponent of Galactic contaminants with plausible (e.g., uniform 
within the accepted range) RV distribution. 

Figure [3]plots the Sextans RV dispersion profile calculated 
using a variety of bin sizes. Aside from setting the global RV 
dispersion value, choice of membership sample has little ef- 
fect on the measured profile. As in K04, the top panels of 
Figure|3]use bins having regularly-spaced outer boundaries at 
5', 10', ...,30'. Our data do not reach the largest radii sam- 
pled by the K04 data, so we can neither confirm nor rule out 
the falling dispersion K04 detect at R > 30'. For R < 30' our 
data indicate a fiat Sextans profile. In contrast, K04 measured 
<jp < 1 km s" 1 for the 7 stars in their sample at R < 5' (see Fig- 
ure 4 of K04). The 40 stars from our sample that occupy this 
region have a p = 8.3+0 g km s" 1 . The profiles shown in Figure 
13^, c) have annuli chosen such that bins contain equal num- 
bers of stars. A "cold" inner profile point begins to emerge 
as the number of bins increases, but at the cost of diminished 
statistical significance. The profiles in (c) use 43 bins (6 stars 
per bin), the smallest number for which we find any annulus 
with (j p < 1 km s" 1 . Monte Carlo simulations indicate that 
one expects to measure such a small dispersion in ~ 1 of 43 
bins even if the true profile is perfectly flat. 5 



We conclude that our data are consistent with a flat Sex- 
tans RV dispersion profile for R < 30'. This result is consis- 
tent with K04 if one omits their innermost point, as K04 do 
when using their observed profile to estimate a Sextans mass 
of M(r < 1 kpc) ~ 10 7 — 10 8 M . We require a more radi- 
ally extended data set in order to further constrain the Sextans 
mass. 

Our data set consists of measurements 
{(X u Y u V u ai),...,(X N} Y N ,V N ,a N )}, where X/ and Y,- 
give the projected position of the i ,h member star with respect 
to the Sextans center, and V, and cr, are the measured RV 
and its formal uncertainty. To identify potential kinematic 
substructure we search for regional anomalies in the RV 
distribution as a function of projected position. From the 
data we estimate the Sextans RV dispersion surface using a 
non-parametric Nadaray a- Watson (1964) estimator: 



x-X; y-Y, 



y-f; 



(1) 



It is perhaps intriguing that our six stars contributing to the cold inner 



where (V) is the sample mean, Ki'^ 1 ,^-) is a smoothing 

kernel and h x , h y specify spatial smoothing bandwidths. This 

estimator provides a weighted local average of (v 2 ) on a grid 
of (x,y). For simplicity we adopt a bivariate Gaussian kernel 

with isotropic smoothing: K(^-, 2zS) oc exp[-| (A '~ A " l) "^ b '~ i ' ) ], 
where h s = h x = h y . Since the data points have a nonuniform 
spatial distribution, we allow the spatial smoothing band- 
width, h s , to vary over the grid surface. At every (x,y), h s 
takes the minimum value for which there exist n data points 
located within 3h s of (x,y). The value chosen for n determines 

point in Figure |3j are localized — all lie within 3' of the Sextans center and 
five occupy a region of diameter 2'. However, Gaussian-random samples of 
294 stars contain, on average, 3 mutually exclusive sets of six neighboring 
stars for which we would measure a p < 1 km s . Including the stars con- 
tributing to the cold point in Figure|3j;, we detect two such sets in our sample. 
As they are consistent with what is expected to arise by chance, we conclude 
that these detections alone do not provide compelling evidence of cold sub- 
structure. 



Kinematic Substructure in Sextans 
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the number of neighboring data points that contribute signifi- 
cantly to the kernel. 
We use Equation ^ to estimate the Sextans RV disper- 
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sion surface, (v 2 ) s (x,y), repeating for neighbor values n = 
10,20, 150. In all cases we evaluate Equation [2 on a 
square grid with x and y values each spanning the range 
{-20',-19',...,+19',+20'}. The top row of panels in Fig- 
ure|4]maps contours of the Sextans RV dispersion surface re- 
sulting from a subset of these estimations with N=294 and 
n = 20,40, 60, 80. 

To assess the significance of apparent surface features we 
adopt as a null hypothesis that the measured RVs are drawn 
everywhere from the same (not necessarily Gaussian) distri- 
bution. To simulate the null hypothesis we generated (for each 
of the three Sextans samples) 1000 artificial data sets in which 
the stars have identical positions to those in the Sextans mem- 
ber sample, but in which (V,-,<7,) pairs are drawn randomly, 
without replacement, directly from the Sextans data. All the 
observed Sextans RVs are present in each artificial data set, 
but scrambled with respect to position. This dissociates any 
existing correlation between position and RV while preserv- 
ing the global RV distribution. 

We apply Equation [2 to each artificial data set and record 
the distributions of artificial surface minima and maxima, 
{v 2 ) Amin and (v 2 } Amax . The significance of a small Sextans 
surface dispersion at grid point (x,y) is then given by the frac- 
tion, p co id(x,y), of artificial surfaces for which (v 2 ) s (x,y) < 
(v 2 ) Amj ; conversely, the significance of a large Sextans dis- 
persion at (x,y) is given by the fraction, pi l0t (x,y), of artificial 
surfaces for which (v 2 ) s (x,y) > (v 2 ) Amax - The bottom two 
rows of panels in Figure |4] map contours of p co id(x,y) and 
Phot(x,y) corresponding to each Sextans RV surface estima- 
tion. 

For nearly all n and all three cases of membership, the Sex- 
tans RV dispersion surface has a minimum value near a point 
~ 15' north of the Sextans center. For the default N=294 
member sample the significance is p co id > 0.9 for estimations 
in which 30 < n < 80, and reaches p co u — 0.965 for n = 40. 
The same northern region remains the coldest in estimations 
using the N=276 and N=303 member samples, despite the fact 
that the former sample removes 2 stars from within 10' of the 
region's center, and the latter contributes 3 high-velocity stars 
to the region. For N=276, the significance is p co id > 0.9 for 
a wide range of n and reaches p co id = 0.948 for n = 60. 6 For 
N=303 the northern feature remains the Sextans surface mini- 
mum, but with only p co id < 0.671. The RV distribution at the 
Sextans center does not provide compelling evidence against 
the null hypothesis. For N=276 and N=294, all surface points 
inside R < 5' have p co id < 0.25; for N=303 the center has 
Pcoid < 0.5, for all n. 

We have also employed the parametric test described by 
Kleyna et al. (2003; K04) which asks whether subregions of 
the position-RV data are more likely described by a single 
Gaussian or by a mixture of two distinct Gaussians. This test 
yields results similar to those from our nonparametric test. 
For the N=294 sample, a model in which RVs from 80% of 
the 35 stars centered on a location 15' NNE of the Sextans 

6 The N=276 Sextans sample contains a second local minimum, 10' south 
of the Sextans center, that is significant with p co M ~ 0.95 in estimations using 
20 < n < 40. Unlike the northern feature, however, this feature disappears 
when considering less restrictive membership cases. 



center (dashed circle in Figure|2j are drawn from a secondary 
Gaussian with disperion 2.7 km s" 1 has likelihood 3.7 x 10 4 
times larger than a global Gaussian model. Simulations using 
artificial data sets indicate this result is sufficient to rule out 
our null hypothesis with significance p = 0.968 (here, p indi- 
cates the fraction of artificial data sets for which no subregion 
yields as large a likelihood ratio). For N=276, the northern 
region is kinematically cold (80% of 35 stars drawn from a 
secondary Gaussian with 2.8 km s" 1 dispersion) with signifi- 
cance p = 0.944. For N=303, the region remains the strongest 
candidate for cold substructure (60% of 35 stars drawn from a 
secondary Gaussian with km s" 1 dispersion), but the signifi- 
cance falls to p = 0.756. No aperture overlapping the Sextans 
center favors a mixed Gaussian model with p > 0.5, for any 
of the three samples and a variety of aperture sizes. 

The striking agreement between parametric and nonpara- 
metric tests prompts two conclusions: 1) we do not confirm 
the presence of a kinematically cold Sextans core, and 2) we 
find stronger evidence (p ~ 0.95 for N=276 and N=294) for 
cold substructure north of center, near the Sextans core ra- 
dius. The latter result is somewhat sensitive to membership 
selection — when the most marginal members are allowed into 
the sample, the significance of the feature drops markedly. 

3. DISCUSSION 

Kleyna et al. have graciously provided us with their RV 
data for 118 Sextans stars (88 members). Velocity measure- 
ments for the 40 stars common to our survey generally agree 
to within formal uncertainties. For R < 5', the data sets con- 
tain 2 stars in common, with excellent RV agreement. What 
accounts for the discrepant results between the two studies? 
Sample size must play a role. K04 propose the existence of a 
kinematically distinct Sextans core based on the cold RV dis- 
persion of the seven innermost stars in their sample. We find 
that if we restrict our analyses to similar numbers of stars, we 
obtain similar evidence. Our RV dispersion profile (Figure 
Qt) identifies a group of six stars near the Sextans center that 
have <7 P = 0.0+ ( 3 ) 4 km s -1 . If we apply the Gaussian mixture 
test to the five of these stars clustered nearest a point 3' north 
of the Sextans center, a scenario in which four are drawn from 
a secondary Gaussian with km s" 1 dispersion is 115 times 
more likely than the alternative in which all five are drawn 
from the main Sextans distribution (for comparison, K04 cal- 
culate a likelihood ratio of 195 for their cold inner sample). 
At the location of these five stars, the squared RV dispersion 
estimated using Equation [2 with n = 5 is negative, indicat- 
ing that the measurement uncertainties dominate the scatter 
among the local velocities. However, for our data neither test 
result is significant at more than the p = 0.5 level. 

The substructure candidate we detect lies near the Sextans 
core radius. If we assume 80% (the fraction indicated by the 
Gaussian mixture test) of the stars in the 15' spanned by this 
feature belong to a distinct "cluster" of stars, and adopt a sur- 
face brightness over this region that is constant and equal to 
one-half Sextans' central surface brightness, we estimate a 
"cluster" luminosity of 3 x 10 4 L Q . This is similar to the lumi- 
nosity Kleyna et al. (2003) ascribe to the cold clump in UMi, 
but an order of magnitude smaller than the luminosity K04 
estimate for a nominal core cluster in Sextans. The timescale 
for the orbital decay of our "cluster" due to dynamical fric- 
tion would therefore be an order of magnitude larger than the 
0.7-1.5 Gyr K04 calculate for a core cluster (assuming the 
same cluster M/L and Sextans halo parameters). This would 
remove the immediate difficulty in explaining why, if the sub- 
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FIG. 4. — Top row: Contours of the Sextans radial velocity dispersion surface estimate, (v 2 ) (x,y), obtained by applying Equation^to the N=294 member 
sample. From left to right, panels indicate results from estimations using nearest-neighbor parameter n = 20,40,60, 80. Black markers indicate the positions of 
member stars used in the estimation. Middle/Bottom row: Contours of statistical significance for regions of apparent cold/hot kinematic substructure. p co u(x,y) 
is the fraction of artificial data sets for which the simulated surface minimum is greater than the estimate from the Sextans data at (x.y). Phot( x ,y) is the fraction 
of artificial data sets for which the simulated surface maximum is less than the estimate from the Sextans data at (x,y). Advancing from outermost to innermost, 
contours enclose regions where p(x,y) > 0.0,0.1,0.2,..., 1.0. 



structure we detect is associated with a bound cluster, it has 
not spiralled to the Sextans center. If such a cluster exists it 
should be detectable with deep imaging. 

Goerdt et al. (2006) provide an explanation for the presence 
of substructure near the core radius of a dSph. They simulate 
the evolution of globular cluster orbits in the Fornax dSph for 
both cored and cusped inner dark matter halos. For a cusped 
halo, dynamical friction drags clusters to the Fornax center 
on timescales less than 5 Gyr. For a cored halo, the infalling 
clusters "stall" indefinitely near the dSph core radius, a con- 
sequence of what Goerdt et al. interpret as a type of orbital 
resonance. Thus the presence of five globular clusters near 



Fornax's core radius strongly favors a cored dark matter halo 
model. In this context it is intriguing that Ursa Minor and 
perhaps Sextans also show signs of kinematically cold sub- 
structure near their core radii. 
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